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Abstract: 

A precise measurement of the inclusive deep-inelastic e + p scattering cross sec- 
tion is reported in the kinematic range 1.5 < Q 2 < 150 GeV 2 and 3 • 10~ 5 < 
x < 0.2. The data were recorded with the HI detector at HERA in 1996 and 
1997, and correspond to an integrated luminosity of 20 pb -1 . The double differ- 
ential cross section, from which the proton structure function F 2 (x, Q 2 ) and the 
longitudinal structure function F L (x,Q 2 ) are extracted, is measured with typ- 
ically 1% statistical and 3% systematic uncertainties. The measured derivative 
(dF 2 (x, Q 2 )/d\n. Q 2 ) x is observed to rise continuously towards small x for fixed 
Q 2 . The cross section data are combined with published HI measurements at 
high Q 2 for a next-to-leading order DGLAP QCD analysis. The HI data deter- 
mine the gluon momentum distribution in the range 3 • 10~ 4 < x < 0.1 to within 
an experimental accuracy of about 3% for Q 2 = 20 GeV 2 . A fit of the HI mea- 
surements and the jip data of the BCDMS collaboration allows the strong cou- 
pling constant a s and the gluon distribution to be simultaneously determined. A 
value of a s (M|) = 0.1150 ± 0.0017(ea;p)lHoo5 (model) is obtained in NLO, 
with an additional theoretical uncertainty of about ±0.005, mainly due to the 
uncertainty of the renormalisation scale. 
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1 Introduction 



Deep-inelastic lepton-nucleon scattering (DIS) has been pivotal in the development of the 
understanding of strong interaction dynamics. Measurements of the inclusive DIS cross 
section have been essential for testing Quantum Chromodynamics (QCD) [Tip. Previous fixed 
target DIS experiments have observed scaling violations, i.e. the variation of the structure 
functions with Q 2 , the squared four- momentum transfer between lepton and nucleon, for 
fixed values of Bjorken-x, which are well described by perturbative QCD. The Q 2 evolution 
of the proton structure function F 2 (x, Q 2 ) is related to the gluon momentum distribution in 
the proton, xg(x, Q 2 ) , and to the strong interaction coupling constant, a s . These can be 
determined with precision deep-inelastic scattering cross section data measured over a wide 
range of Bjorken-s and Q 2 . 

The first measurements of F 2 at low x ~ 1CT 3 and Q 2 ~ 20 GeV 2 at HERA revealed 
a steep rise of F 2 (x, Q 2 ) towards low x for fixed Q 2 [@,|3|]. Strong scaling violations are 
observed at low x and are attributed to a high gluon density in the proton. The validity of 
the DGLAP evolution equations [TJfJ, which neglect higher-order terms [|5|, |5p proportional 
to a s ■ ln(l/x), is questionable in the low x range and therefore has to be tested against 
data. At extremely low x, non-linear gluon interaction effects have been considered in order 
to damp the rise of the cross section in accordance with unitarity requirements [|7p. The 
study of quark-gluon interaction dynamics at high parton densities therefore continues to 
be a challenging subject. Knowledge of the parton densities at low x is also necessary for 
interpreting measurements at hadron colliders and of cosmic neutrino interactions. 

This paper presents new cross section measurements for the neutral current process 
e+p -> e + X in the kinematic region 1.5 < Q 2 < 150 GeV 2 and 3 ■ 10~ 5 < x < 0.2, ob- 
tained from data taken in the years 1996 and 1997 with positrons of energy E e = 27.6 GeV 
and protons of energy E p = 820 GeV, corresponding to a centre of mass energy */s = 300.9 
GeV. Cross section measurements at low x and medium Q 2 , based on the 1994 HERA data, 
were previously published by the HI collaboration [|8p and by the ZEUS collaboration [|9|]. 
The present measurement uses upgraded detectors to measure and identify the scattered 
positron, including new precision tracking for low Q 2 scattering. It also benefits from in- 
creased luminosity from HERA which enables an accuracy of typically 3% to be reached for 
the DIS cross section. Thus it considerably improves the former structure function measure- 
ments of the HI collaboration [||,[T^,|TT]] at Q 2 < 150 GeV 2 . The kinematic range is extended 
to larger x, yielding an overlap of HI data with measurements from fixed target muon-proton 
scattering experiments for the first time. The paper includes a measurement of the derivative 
(dF 2 /d\n Q 2 ) x , which serves as a sensitive test of the dynamics of strong interactions. 

The longitudinal structure function F L (x, Q 2 ) is obtained with improved precision in an 
extended range of inelasticity y and Q 2 , as compared to its first determination at low x JTT]] . 
A new method is introduced to extend the extraction of Fl to Q 2 values below 10 GeV 2 which 
uses the derivative (da r /d In y)q2 of the reduced cross section oy. 

A next-to-leading order (NLO) DGLAP QCD analysis is performed using inclusive lepton- 
proton scattering data by introducing a new flavour decomposition of the structure function 
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F 2 . Hence it is independent of nuclear binding effects in the deuteron or heavier nuclei. The 
QCD analysis of the present low x data and of the recently published high Q 2 HI data flT2l ] 
determines the gluon distribution xg(x, Q 2 ) at low x. The combination of the low x HI data 
with large x data from the BCDMS experiment JT3| ] enables an accurate, simultaneous de- 
termination of both xg(x, Q 2 ) and a s (M|) . The present analyses use all of the available 
information regarding the experimental uncertainties of the data sets considered and explore 
the QCD model and fit parameter variations in a systematic way. 

The paper is organised as follows. Section 2 defines the inclusive cross section and the 
methods used to reconstruct the event kinematics. The detector, the event selection and the 
simulation are described in section 3. Section 4 presents the alignment and calibration meth- 
ods, and summarises the cross section measurement. Section 5 presents the measurement of 
(da r /dhiy)Q2 and the determination of the longitudinal structure function F L . The results 
for the proton structure function F 2 and its derivative (dF 2 /d In Q 2 ) x are given in section 6. 
The QCD interpretation of the data is discussed in section 7, which refers to an appendix 
presenting details of the analysis. The paper is summarised in section 8. 



2 Cross Section and Kinematic Reconstruction 

The inclusive DIS cross section of the measured reaction e + p — > e + X depends on two in- 
dependent kinematic variables, chosen to be x and Q 2 , and on the centre of mass energy 
squared s, with the inelasticity variable y = Q 2 / sx. In the one-photon exchange approxima- 
tion the neutral current double differential cross section, d 2 a/dxdQ 2 , is given by 



d 2 a 2na 2 Y, 



dxdQ 2 Q 4 x 
where the reduced cross section is defined as 



°r (1) 



a r = F 2 (x,Q 2 )-^-F L (x,Q 2 ) (2) 

and Y + = 1 + (1 — y) 2 . Due to the positivity of the cross sections for longitudinally and trans- 
versely polarised photons scattering off protons, the two proton structure functions F 2 and 
F L obey the relation < F L < F 2 . Thus the contribution of the longitudinal structure func- 
tion F L to the cross section can be sizeable only at large values of the inelasticity y, and in 
most of the kinematic range the relation a r « F 2 holds to a very good approximation. 

The HERA collider experiments allow DIS kinematics to be reconstructed using the scat- 
tered positron, the hadronic final state, or a combination of the two. This is important for 
maximum coverage of the kinematic range and the control of systematic uncertainties. 

In the "electron method" the event kinematics are determined using the measured energy 
of the scattered positron E' e and its polar angle 9 e according to the relations 

rpf J?' 2 air, 2 A 

Ve = 1 - £ sin 2 (^/2), Ql = \ Ue . (3) 



The coordinate system of HI is defined such that the positive z axis is in the direction of the 
incident proton beam. Polar angles 9 are defined with respect to the proton beam direction. 
While the electron method is accurate at large y, corresponding to low E' e , the resolution 
rapidly degrades with l/y as E' e approaches the positron beam energy E e . 

The inelasticity y can also be determined as [ P"4"[ ] 

_ Xi(Ei-p Zii ) = _£_ 
Vh 2E e ~2E e ' K) 

where Ei and p z>i are the energy and longitudinal momentum component of a particle i in 
the hadronic final state, the masses being neglected. In this analysis the kinematics are also 
reconstructed with the "S method" using the variables JT^ ] 

£ _ 2 _ E' e 2 sin 2 fl e 

E + ££(1 - cos6> e ) 1-2/s 

For all reconstruction methods, Bjorken-x is calculated as i = Q 2 /sy. Due to energy- 
momentum conservation the variable 

E-p z = Z + E' e (l-cos9 e ) (6) 

is approximately equal to 2E e . The hadronic variables y h and y s are related according to 

ys = — (7) 

i + y h -y e 

and can be well measured down to low y ~ 0.004. The variable ys is less sensitive to initial 
state radiation than y h since the initial energy E e in the denominator in equation |] can be 
calculated using the total energy reconstructed in the detector which leads to equation |5|. 

The hadronic scattering angle is defined as 

tan| = -^, (8) 

where P t:h is the total transverse momentum of the hadronic final state particles. In the naive 
quark parton model, 9 h defines the direction of the struck quark related to 9 e as 

tan^ = ^.tan^. (9) 
2 1-y 2 

For y > 0.5 the positron scattering angle 9 e is smaller than 9^. This relation, together with 
the definition of y e (equation 0), determines the scattered positron energy from 9 e and 9h in 
the "double angle method" 



6 



3 Experimental Procedure 



3.1 HI Detector 

The HI detector [ fT7| ] combines tracking in a solenoidal magnetic field of 1.15 T with nearly 
hermetic calorimetry to investigate high-energy ep interactions at HERA. The low Q 2 cross 
section measurement relies mainly on the central and backward tracking systems, the back- 
ward calorimeter (SPACAL) and the Liquid Argon (LAr) calorimeter. These components are 
briefly described below. 

The energy of the positron, when scattered into the backward region of the HI detector 
(153° < 9 e < 177°), is measured in the SPACAL, a lead-fibre calorimete r [jTg,|Tgfl . The 
SPACAL has an electromagnetic section with an energy resolution of 7%/ JE/GeV © 1%, 
which together with a hadronic section represents a total of two interaction lengths. The 
SPACAL time resolution of less than 1 ns allows proton beam induced background to be 
vetoed. The calorimeter has a high transverse granularity which provides a determination of 
the transverse coordinates of electromagnetic clusters with an accuracy of a few millimeters 
and positron identification capability. Identification of the scattered positron is improved 
and the polar angle measured with a backward drift chamber (BDC), situated in front of 
the SPACAL, and a new backward silicon strip detector (BST) [|20|]. The BST consists of 
four detector planes, arranged perpendicularly to the beam axis which are equipped with 16 
wedge shaped, double metal silicon strip detectors. The BST measures the polar angle of 
tracks with an internal resolution of about 0.2 mrad from radial coordinates between 5.9 cm 
and 12.0 cm. 

The hadronic final state is reconstructed with the Liquid Argon calorimeter [^fl], the 
tracking detectors and the SPACAL. The LAr calorimeter is built of eight wheels of modules 
with an octant structure. The total depth of the calorimeter varies between 4 and 8.5 interac- 
tion lengths depending on the polar angle. In test beam measurements pion induced hadronic 



energies were reconstructed with a resolution of about 50%/y E /GeV (B2.0% after software 
energy weighting 

The interaction vertex is determined with the central tracking detector consisting of two 
concentric jet drift chambers (CJC) and two concentric z drift chambers (CIZ and COZ). The 
vertex determination is complemented by the inner proportional chamber CIP, for 167° < 
6 e < 171°, and by the silicon tracker BST, for 171° < 6 e < 176.5°. 

The luminosity is determined using the small-angle bremsstrahlung process ep — > ep 7. 
The final state photon and the positron, scattered at very low Q 2 , can be detected in calorime- 
ters ("electron and photon taggers") which are situated close to the beam pipe at distances of 
33 m and 103 m from the interaction point in the positron beam direction. The luminosity is 



measured with a precision of 1.5% using the method outlined in [23]. 



3.2 Data Samples and Interaction Triggers 



The analysis comprises two different DIS event samples: 
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• Sample A - data taken in the years 1996 and 1997 with luminosities of 4.5 pb 1 and 
13.4 pb -1 , respectively. These two data sets are combined to provide the cross section 
measurement for Q 2 values from 15 GeV 2 to 150 GeV 2 and for Q 2 = 12 GeV 2 at y > 
0.17. 

• Sample B - data taken in the autumn of 1997 during a two week period dedicated to the 
lower Q 2 region. The data from this special run with a luminosity of 1 .8 pb -1 are used 
in the Q 2 range from 1.5 GeV 2 to 8.5 GeV 2 and for Q 2 =12 GeV 2 at low y < 0.17. 

DIS events at low Q 2 are characterised by a positron scattered into the backward part of the 
HI apparatus. The event trigger for sample A requires local energy sums in the SPACAL 
calorimeter to be above an energy threshold of 6 GeV. This threshold is lowered to 5 GeV in 
sample B. 

Both data samples are contaminated by photoproduction events in which the scattered 
positron escapes through the beam pipe and a particle in the hadronic final state mimics the 
signature of a scattered positron. For a fraction of these background events the scattered 
positron is detected in the electron tagging calorimeter. This background is significant only 
at low energies E' e , corresponding to values of y > 0.6. 

The region of high y > 0.75 is accessed with a dedicated trigger which requires a com- 
pact energy deposition (cluster) of more than 2 GeV of energy in the SPACAL, and a vertex 
signature in the proportional chamber system. The data accumulated with this trigger corre- 
spond to a luminosity of 2.8 pb -1 in 1996 and 3.4 pb 1 in 1997. 

The SPACAL energy triggers are monitored with independent track triggers and found to 
be fully efficient for energies of about 1 GeV above threshold. The high y trigger efficiency is 
determined to be 97% using independent calorimeter triggers. The online data reconstruction 
leads to a maximum loss of 0.5% of DIS events. This loss is estimated from monitor data 
and corrected for. 



3.3 Event Selection 

The scattered positron is identified with the cluster of maximum transverse momentum p t in 
the SPACAL calorimeter, for which requirements on the cluster shape are satisfied. Electro- 
magnetic energy deposition leads to clusters of smaller transverse extension than hadronic 
energy deposition. The transverse energy distribution of positron showers is determined ex- 
perimentally using QED Compton events, and from radiative DIS events in which a photon is 
radiated from the incoming positron and detected in the photon tagger. The positron cluster 



radius Q19|,124|] can thus be measured at all energies considered in the analysis, and a cluster 
radius cut of 4 cm is chosen. This cut removes a sizeable fraction of the photoproduction 
background while retaining more than 99% of the DIS signal. 

A positron candidate cluster is required to be associated with a track segment in the BDC. 
The efficiency of the BDC is measured to be 98% on average with small radius dependent 
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variations. A signal is also required in either the CIP or the BST where geometrically avail- 
able. The efficiency of the CIP is about 98%. The track efficiency of the BST in the 1997 
special run period is found to be about 93%. Efficiencies at low energies are determined 
using QED Compton scattering events and radiative DIS events. 

Reconstruction of the interaction vertex is necessary to determine the event kinematics 
and to suppress beam background events. In the intermediate y region, hadrons measured in 
the central tracking chambers allow vertex reconstruction with an efficiency exceeding 98%. 
However, at very low y and also at very high y, where no hadron may be measured in the 
central tracking chambers, the vertex can be defined by the scattered positron if it falls within 
the acceptance of the CIP or the BST. 



SPACAL energy 


> 6.9 GeV 


radius of cluster 


< 4 cm 


fraction of energy in the hadronic section 


< \5%ofE' e 


cluster-BDC track distance 


< 1.5 cm 


cluster-BST track distance 


< 1.0 cm 


z vertex position 


\z\ < 30 cm 


E-Vz 


> 35 GeV 



Table 1: Basic DIS event selection criteria. The radius of the cluster /[I^,|2^] defines the 
lateral shape of the energy deposition. The high y region, corresponding to energies between 
3 GeV and 6.9 GeV, is accessed with additional cuts, as discussed in Section PO. 



K 



Longitudinal momentum conservation in neutral current DIS events gives the constraint 
that E — p z , summed over the final state particles is about equal to 2E e (equation In 
radiative events a photon may carry a significant fraction of the E — p z sum. Such events are 
thus suppressed by requiring E — p z > 35 GeV. 

The criteria applied to select DIS events are summarised in table [j]. 



3.4 Simulation 



For the calculation of the detector acceptance and efficiency control about 10 7 inelastic events 



are simulated. Deep-inelastic events are generated using the DJANGO [ |25| ] event generator. 



This program is based on HERACLES [26] for the electroweak interaction and on the Lund 



Monte Carlo generator program ARIADNE Q27Q, which includes the generation of diffrac- 
tive events. This generator, when tuned to HERA data [[281], presently gives the most reliable 



description of the final state properties [|29|]. To describe higher order QCD radiation pro- 
cesses the ARIADNE program uses the Colour Dipole Model (CDM) pD]]. For hadron 
fragmentation the JETSET program is used [[3lj]. Comparisons are done using the generator 
LEPTO [g2p. QED Compton events are generated using the program COMPTON [T53]]. DIS 
events are also generated with the HERWIG event generator [|34]], which includes resonant 
final state production. This is important for the description of rejected events at very low 
y p5|]. Photoproduction background is generated with the PHOJET [P6fi program using the 
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parameterisation of CKMT [ p7| ] to determine the virtual photon-proton interaction cross sec- 
tion. The normalisation of the PHOJET event sample is adjusted to the data measured with 
the electron tagging calorimeter. It is found to agree within 20% with the calculation of the 
cross section using the Weizsacker- Williams approximation. Using the leading logarithmic 
approximation the effect of photon radiation is estimated to be negligible. 



The detector response is simulated in detail with a program based on GEANT3 H39fl. 
The Monte Carlo events are subject to the same reconstruction and analysis chain as the 
real data. In the comparisons shown here, the simulated distributions are normalised to the 
measured luminosity. In the event generation the DIS cross section is calculated with the 
parton distributions of p0| ] and with the longitudinal structure function — 0. A QCD fit 
to all the data is used to reweight the simulated cross section. 



4 Measurement of the Cross Section 



This section presents briefly the methods and results of the measurement of the deep-inelastic 
scattering cross section. Further details of the analysis are described in [ ^4||35| , pd1 , |42| ]. 



4.1 Detector Alignment 

The coordinate system of the HI detector is defined by the central tracking chambers which 
determine the spatial coordinates of the interaction vertex. The variation of the vertex po- 
sition in the x, y plane along the beam direction is used to determine the inclination of the 
beam axis with respect to the z axis. 

The central tracking chambers are aligned with respect to each other using cosmic muon 
tracks. The alignment of the BDC and of the SPACAL with respect to the central tracker 
is done by studying the difference of the polar angles measured by these detectors as a 
function of the azimuthal angle. This results in 1 to 2 mm adjustments of the nominal detector 
positions. QED Compton events, which have the signature of back-to-back positron and 



photon clusters, provide a cross check for the alignment Q24J ] of the SPACAL in the transverse 
plane to an accuracy of 0.2 mm. After internal adjustment of the strip detector planes, the 
spatial position of the BST is determined using the event vertex z coordinate measured with 
the central tracker. 

In the BST angular acceptance range the polar angle is measured both by the BDC using 
the interaction vertex, and by the BST track segment. This allows the alignment procedure 



to be checked to within an accuracy of A9 of 0.1 mrad Q35[] . From the residual dependence 
of A9 on the azimuthal angle and from the uncertainties of the alignment procedure, a mea- 
surement error of 0.3 mrad is estimated for the angle of the scattered positron. 
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4.2 Calibration of the Energy Measurements 



The energy of the scattered positron is measured in the SPACAL, which has a transverse cell 
size of 4 x 4 cm 2 and a Moliere radius of 2.5 cm. In a first step the responses of the SPACAL 
cells are equalised using cosmic muons. The energy scale of each cell is determined with DIS 
events using the double angle method, which allows the energy of the scattered positron to 
be expressed as a function of 9 e and 8h- This method is applied to the data of the various run 
periods, and also to the simulated events. Agreement of the energy scales is found at the level 
of 0.2%. The calibration procedure leads to a systematic error of the E' e scale of 0.3% for 
most of the SPACAL area and energies E' e above 20 GeV. The energy calibration at lower 
energies is performed using QED Compton events. This leads to an estimated maximum 
energy scale uncertainty of 2.7% at 3 GeV. This uncertainty is observed to approach linearly 
the 0.3% level at maximum energies E' e ~ E e . The SPACAL response at the lowest energies 
is cross checked by studying the n° mass reconstructed from pairs of photons in the energy 
range of 0.8 — 4 GeV [^3],|3p. Material in front of the calorimeter leads to showering and 
energy losses, which are corrected for using the backward drift chamber as a preshower 
detector [|TJ. Figures [T]a,b and ^a,b show the energy and polar angle distributions for data 
samples A and B. These are described by the simulation of DIS and photoproduction events. 

The cross section measurement at low y relies on the measurement of y using the hadronic 
final state (equation |]). The determination of y is optimised by combining calorimeter energy 
deposits with low momentum tracks. The sum over energy clusters in the calorimeters can be 
strongly affected by electronic noise, in particular for low y < 0.03. Thus an additional noise 
suppression is performed which excludes isolated depositions of energy less than 400 MeV 
(800 MeV) in the central (forward) region of the LAr calorimeter from the analysis of both 
data and simulated events. This leads to a small signal loss but improves the y resolution at 
low y. The uncertainty of this subtraction procedure is estimated to correspond to a quarter 
of the suppressed noise contribution to the reconstructed hadron energy. 

The calibration of the hadronic energy measurement uses the p t balance between the 
scattered positron and the hadronic final state. The energy scales for the electromagnetic and 
hadronic sections of wheels and octants [ |35| ] of the LAr calorimeter result from a Lagrangian 
multiplier technique which simultaneously determines all 128 calibration constants. A sys- 
tematic uncertainty of 2% on the hadronic energy scale in the LAr calorimeter is determined, 
based on a binwise comparison of different calibration methods. The results of the energy 



calibration procedure are consistent with those recently presented [12]. 



Figures [T|c and ||c show the experimental y^ distributions over three orders of magnitude 
in y. In the simulation, the y^ reconstruction is found to be accurate to within a few % 
over a wide range of y h and down to the small Q 2 region. The worse y resolution at low 
y is accounted for with an increased bin size, allowing the y range to be extended down to 
y ~ 0.004. 

The response of the SPACAL to hadronic energy flow is calibrated using longitudinal 



momentum conservation in the DIS events [ |24| ] to within an uncertainty of 5%. This scale 



uncertainty affects the final cross section data at large y through the E — p z momentum 
balance requirement. 
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4.3 Measurement at Large y 



For the measurement of the longitudinal structure function it is essential to reach the max- 
imum possible values of y (see equation Q). This requires an efficient rejection of photo- 
production background events in which low energy deposits in the SPACAL can mimic the 
signature of a deep-inelastically scattered positron. 

At Q 2 below 10 GeV 2 , the range y < 0.75 is accessed by requiring a track signal in 
the BST. This requirement removes a sizeable fraction of the background where a cluster in 
the SPACAL is due to photons from n° — > 77 decays. The remaining background is due to 
photon conversion and showering in the passive material, possible overlap of ir° decays with 
charged tracks, and misidentified charged pions. This background is subtracted bin by bin 
using the PHOJET event simulation. Figures |]a and illustrate the range of polar angle and 
energy for high y events with a track in the BST. The photoproduction background can be 
estimated experimentally using a data sample of events which have a positron detected in the 
electron tagger. Figure shows the energy spectrum of SPACAL clusters for those events 
which satisfy the DIS selection criteria, apart from the E —p z requirement. This distribution 
is well described by the simulation. 

At Q 2 above 10 GeV 2 , for y < 0.75, the photoproduction background is subtracted 
using the PHOJET event simulation. For y above 0.75, however, experimental information is 
used by employing the charge assignment of central tracks associated with SPACAL energy 
clusters. This allows the energy range to be extended down to E' e = 3 GeV, corresponding 
toy < 0.89. For 12 < Q 2 < 35 GeV 2 tracks reconstructed in the CJC can be linked to 
low energy SPACAL clusters with an efficiency of 95% (93)% in 1996 (1997). For such 
tracks with energies up to 15 GeV the charge is determined with an efficiency of 99.5% p3fl. 



The sample of candidates with negative charge is taken to represent the background in the 
positron data sample. 

The statistical subtraction procedure requires the study of any process which may cause 
a charge asymmetry. This asymmetry can be measured using tagged photoproduction events 
which fulfill the DIS event selection criteria. A small charge asymmetry (N + — AL) / (N + + 
NJ) is found with an average of —4.8% with a statistical accuracy of 1.9%, for 0.65 < y < 
0.89. Here N + (N_) is the number of events with positive (negative) charge of the track 
associated with the SPACAL cluster. Comparing the energy distribution for a sample of neg- 
ative tracks in e + p scattering, taken in 1996/1997, with that for a sample of positive tracks 
in e~p scattering, taken in 1999, a consistent asymmetry of (—3.5 ± 2.5)% is measured. 
Background simulation studies and measurements of the ionisation losses in the CJC reveal 
that this small asymmetry is due to the antiproton interaction cross section exceeding that 
for proton interactions at low energies [Hfl. Annihilation leads to larger energy deposits in 



SPACAL than proton interactions which introduces an asymmetry for low energies above 
a given threshold. This charge asymmetry is taken into account in the measurement of the 



positron DIS cross section at high y. Selected control distributions [ [24| , |42| ] are shown in fig- 
ure ^ illustrating the good understanding of this kinematic region down to scattered positron 
energies of E' e = 3 GeV. 
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4.4 Results 



An iterative bin wise correction procedure is adopted for the extraction of the double differ- 
ential cross section oy. This procedure requires that the bin sizes are adapted to the resolution 
in the measurement of the kinematic variables. The data and the simulated events are binned 
in a grid in x with five bins per decade and in Q 2 with eight bins per decade, as illustrated in 
figure |5|. At low y the resolution of the measurement degrades and the bin size is widened. 
For y > 0.6 the data are divided in bins of y, and the Q 2 division is kept. In this region 
the cross section may receive a large negative contribution from F L proportional to y 2 and 
therefore a fine binning in y is desirable. Bins are accepted if the purity and stability are 
bigger than 30% with typical values being 60%. Here the purity (stability) is defined as the 
number of simulated events which originate from a bin and which are reconstructed in it, 
divided by the number of reconstructed (generated) events in that bin. 

The longitudinal momentum conservation constraint, E—p z > 35 GeV, limits the amount 
of radiative corrections to at most 5% at high y. The program HERACLES [E5p, which is 



used in the DIS event simulation, accounts for first order radiative corrections to the one 
photon exchange approximation. The radiative corrections are extracted using a high statis- 
tics calculation within the HERACLES Monte Carlo program, and compared with the results 
of the numerical program HECTOR [|33|], which includes higher order and hadronic correc- 
tions. The corrections are found to agree within the statistical accuracy of the radiative event 
simulation of 0.5%. 



The results of the measurement are summarised in tables PHT2|. At y > 0.17 the kinemat- 
ics are reconstructed using the quantities Q 2 and y e . At y < 0.17, where the resolution of 
y e degrades, the variables Q| and y% are used. The error calculation for the measurement is 
discussed below. The full error correlation matrix can be obtained from the HI Collabora- 
tion ffl 



The cross section measurement is shown in figure |as a function of x for different Q 2 . 
Due to the extension of the measurement towards low y, the HI data overlap with data of 
fixed target fxp experiments. The HI measurement agrees well with the fixed target data 
within the uncertainty of about 7 %. The cross section rises towards low x. This rise is 
observed to be damped at the smallest values of x, which is attributed to the longitudinal 
structure function, see section |5]. The cross section can be well described by a QCD fit to the 
data as discussed in section [7]. 



4.5 Systematic Errors 

The statistics of the data presented here exceed 10 4 events in most of the bins. The precision 
for this measurement is dominated by systematic uncertainties of typically 3%. These are 
classified into a global normalisation uncertainty, kinematically correlated errors (S cor ), the 
statistical errors of the data (5 sta ) and uncorrected errors (5 unc ). The uncorrelated errors 
contain the statistical uncertainty of the simulation and further errors due to local systematic 
effects. 
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Table || lists those errors which result in a possible global change of all data points. The 
resulting total normalisation uncertainty of the data is 1.7%. It is dominated by the error on 
the luminosity measurement. 



source 


cross section error [%] 


online event selection 


0.5 


BDC efficiency 


0.5 


trigger efficiency 


0.5 


luminosity measurement 


1.5 



Table 2: Sources and sizes of normalisation errors. 



Energy calibration and alignment uncertainties cause systematic errors which depend on 
the kinematics and introduce correlations between the measured data points. These errors 
are determined using the simulation program and verified by an analytic calculation. They 
are found to be symmetric to good approximation. The uncertainty of the photoproduction 
background simulation is estimated to be 20% . The correlated error sources are listed in 
table |. 



source 


size of uncertainty 


typical cross section error [%] 


scattered positron energy scale 

scattered positron angle 
hadronic energy scale in LAr 
LAr noise 

photoproduction background 


0.3%at££~27.5GeV 
2.7% atE' e = 3GeV 
0.3 mrad 

2% 

25% of noise 
20% of background 


1 

2 

0.5 

2 

max of 5 at lowest y 
3 at large y 



Table 3: Sources of correlated systematic errors and their typical effect on the cross section 
measurement accuracy. 



As a cross check of the positron identification, the scattered positron is also considered to 
be the cluster of maximum energy. When this alternative positron identification criterion is 
used, the cross section changes by less than 1%. This is accounted for in the systematic error. 
Detailed studies using different event generators with differing simulations of the hadronic 
final state verify the reliability of the positron identification procedure within the quoted 
systematic uncertainty 



Uncertainties due to radiative corrections, positron identification and final state simula- 
tion details are treated as uncorrelated systematic errors. The errors introduced by the track 
based background subtraction procedure in the high y data analysis (section ft~3| ) are sum- 
marised in table |j. 

For each kinematic bin the resulting cross section errors are given in the cross section 
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tables 9- 12 



source 


error of cross section [%] 


positron identification 


1 


track charge determination 


0.5 


charge asymmetry 


1 


CJC-SPACAL track link efficiency 


2 


hadronic track requirement in CJC 


1 


high y trigger efficiency 


1 


radiative corrections 


1 



Table 4: Additional sources of systematic errors and effect on the cross section for the ana- 
lysis, in the range 12 < Q 2 < 25 GeV 2 , in which central jet chamber tracks are used to 
measure the charge of the positron candidate. 



5 Longitudinal Structure Function F L (x,Q 2 ) 



The extraction of the longitudinal structure function is based on the reduced double differen- 
tial cross section (equation §), which depends on two proton structure functions, F 2 (x, Q 2 ) and 
Fl(x, Q 2 ) . The contribution of Fl is enhanced with y 2 , and the reduced cross section cr r 
tends to F 2 — Fl for y — > 1. In the quark-parton model, the longitudinal structure function is 
zero for spin 1/2 quarks p7|]. In QCD, parton radiation processes pS| ] lead to non-zero val- 



ues of F L . Thus F L contains information about the gluon distribution and about the strong 
interaction dynamics which is complementary to that obtained from the analysis of the scal- 
ing violations in F 2 (x, Q 2 ) . At low Q 2 the longitudinal structure function is expected to be 
particularly sensitive to higher-order corrections to DGLAP QCD 



The longitudinal structure function can be extracted from the inclusive cross section only 
in the region of large y. An important advantage of HERA, compared to fixed target DIS 
lepton-nucleon experiments, is the wide range of y values covered. This allows the behaviour 
of F 2 at low y to be determined reliably and to be extrapolated into the region of high y. Two 
methods are used here to perform an extraction of the longitudinal structure function. For 
larger Q 2 , a NLO DGLAP QCD fit is used to extrapolate F 2 into the high y region, and F L is 
obtained with the "extrapolation method" introduced previously [JTTJ] . This fit uses only HI 
data in the restricted kinematic range y < 0.35 and Q 2 > 3.5 GeV 2 . Details of this and 
other QCD fits to HI and fixed target data are presented in section [7]. At low Q 2 < 10 GeV 2 , 
the behaviour of F 2 as a function of In y is used in a new extraction method [|42"|], termed the 
"derivative method" as it is based on the cross section derivative (da r /d\ny)Q2 . 



5.1 Cross Section Derivative (dcr r / d In y)cp 

The derivative of the reduced cross section, taken at fixed Q 2 , is given by 



dlnyJ Q2 \d\ny) Q2 Y 2 d\ny Y. 
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For y — ► 1 the cross section derivative tends to the limit (dF 2 /dlnyjqz — 2 ■ F L , neglecting 
the contribution from the derivative of Fl. At largest y the Fl contribution dominates the 
derivative of the reduced cross section a r . This is in contrast to the influence of F L on a r 
which is dominated by the contribution of F 2 for all y. A further advantage of the derivative 
method is that it can be applied down to very low Q 2 ~ 1 GeV 2 where a QCD description of 
F 2 (x, Q 2 ) is complicated due to higher order and possible non-perturbative corrections. 

To obtain an accurate determination of (dcr r /din y) Q 2 the data are rebinned in Q 2 by 
combining data of two adjacent Q 2 intervals. Differences Aa r are calculated between cross 
section points adjacent in y at fixed Q 2 . A bin-centre correction is applied to obtain the 
derivative at each y point, which is chosen to be the average of the two y values of the cross 
section measurements used to calculate the derivative. A full error analysis is performed in 
order to account for the correlations of errors, which partially cancel. The two adjacent data 
points of the derivative are anti-correlated since they use the same cross section measure- 
ment with different sign. The cross section derivatives are shown in figure [7], and the values 
are given in tables [T3] and [T3]. The measured derivatives are well described by the QCD 



calculation (section |7T2"|). 

For low Q 2 and y < 0.3 the derivative is observed to be a linear function of lay. The 
structure function F 2 (x, Q 2 ) , at fixed Q 2 , behaves like x~ x oc y x . At low Q 2 the exponent A 
is observed to be small and the derivative is thus expected to be approximately linear in In y. 
This approximation is used to determine the longitudinal structure function at low Q 2 . For 
larger Q 2 the exponent A rises and a curvature is expected as can be seen in figure M. 



5.2 Determination of Fl 

For the determination of F L for Q 2 < 10 GeV 2 straight line fits are performed in In y to the 
derivative (da r /d\ny)Q2for y < 0.2. These straight lines describe the data well (figure ^) 
and are extrapolated to estimate the contribution of dF 2 /d\ny at high y. The uncertainty 
on this extrapolation is included in the systematic errors of the Fl determination, taking 
into account the correlations of errors at low y with those at high y. The extrapolations are 
compared with the values obtained from the QCD analysis and very good agreement has 
been found. The contribution of dF L /d\ny to the derivative (equation |T0|) is neglected. 
The uncertainty of this approximation is estimated using the derivative of Fl as calculated 
in QCD. It is taken as an additional uncertainty on the measurement, and amounts to about a 
quarter of the systematic uncertainty of Fl ■ 

For Q 2 > 10 GeV 2 , the NLO QCD fit of the HI data for y < 0.35 is used to estimate F 2 in 
the high y (i.e. small x) region. In this Q 2 range, the extrapolation method is more accurate 
than the derivative method. The derivative method is statistically limited at maximum y since 
this region is accessed only with a combined SPACAL and track trigger. Figure [8| compares 
the fit with the measured cross section for those five Q 2 bins above 10 GeV 2 which access 
the high y region. The difference between the measured a r and the extrapolated F 2 is used to 
determine F L (x, Q 2 ) as described in [|1 iy. Systematic errors, which are common to the lower 
y and the large y region, are considered in the fit as described in [|52|] . 
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The Fl values obtained are presented in table [15]. The uncertainties on the longitudi- 
nal structure function include several sources: the statistical errors, uncorrelated systematic 
errors and correlated systematic errors, resulting e.g. from the y dependent amount of sub- 
tracted photoproduction background. In addition, errors are associated with the assumptions 
inherent to the extraction methods. For the derivative method these errors are dominated by 
the uncertainty of the straight line fit, and for the extrapolation method by the variation of 
the smallest Q 2 of data used in the QCD fit, see section As can be seen in table [13], 
these errors are smaller than the experimental systematic errors. In the region of overlap, for 
Q 2 between 4 GeV 2 and 15 GeV 2 , the derivative method and the extrapolation method give 
consistent results. Further details of this analysis are described in [|42|]. 

The values for Fl(x, Q 2 ) , given in table O, are consistent with the previous determina- 
tion of -Fiby the HI collaboration [JTTjl, but they are more precise and obtained in a wider 
kinematic range. The HI data extend the knowledge of the longitudinal structure function 
beyond that obtained from fixed target lepton-proton scattering experiments into the region 
of much lower x, see figure |^. The increase of Fl(x, Q 2 ) towards low x is consistent with 
the NLO QCD calculation (section reflecting the rise of the gluon momentum distri- 



bution towards low x. A measurement of the x dependence of Fl(x, Q 2 ) can be performed 
independently of assumptions about the behaviour of F 2 with a variation of the proton beam 



energy at HERA [53]. 



6 Structure Function F 2 and Derivative (dF 2 /d In Q 2 ) 



X 



The proton structure function F 2 (x, Q 2 ) is obtained from the measured reduced cross section 
using equation rewritten as 



,2 r, \ -1 



The ratio R = F L /(F 2 - F L ) is determined using the standard DGLAP QCD fit to the 
HI data (section [7T2"| ), calculating Fl to order a 2 . In order to reduce the dependence of the 
measurement on F L , the F 2 extraction is limited to the range y < 0.6. The results for 
F 2 (x, Q 2 ) and the calculated values of F L (x, Q 2 ) are given in tables |9|-|r2|. This measurement 



is consistent with, and improves upon, the previous results [g], which were obtained with a 
different backward apparatus. 

In figure [K] the measurement of the structure function F 2 (x, Q 2 ) at low x is shown as a 
function of Q 2 . The data are well described by the NLO QCD fit, as is discussed in detail 
in section [7T2[ The In Q 2 dependence of F 2 is observed to be non-linear. It can be well 
described by a quadratic expression 

P 2 (x,Q 2 ) = a(x) + b(x) In Q 2 + c(x) (In Q 2 ) 2 , (12) 

which nearly coincides with the QCD fit in the kinematic range of this measurement. 
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The DGLAP evolution equations are governed by the derivative (dF 2 /d In Q 2 ) x taken at 
fixed x. Measurement of this derivative has long been recognised as a powerful constraint on 
xg and a s []53J]. In the low x DIS region its behaviour is a direct reflection of the behaviour 
of the gluon density []55[], This quantity has also been studied in [ J56] ] in view of possible 
non-linear gluon interaction effects O]. A study of the derivative (dF 2 /d\n Q 2 ) x at low x 
was presented previously by the ZEUS Collaboration [ j57] ] where -F 2 (a;, Q 2 ) was assumed to 
depend linearly on \nQ 2 . This approximation is inconsistent with the data presented here, 
see figure [HJ 

Using the procedure adopted for the derivatives (da r /d\ny)Q2 (section |5T| ) the local 
derivatives (dF 2 / dhiQ 2 ) x are measured. The results are shown in figure [IT] for different 
a; as a function of Q 2 , and the values are quoted in tables [16|,[F7] and [18]. For each bin of x 
these derivatives can be described by the function b(x) + 2 • c(x) In Q 2 (solid lines). Small 
deviations of (dF 2 /d\nQ 2 ) x from the linearity in In Q 2 occur in NLO QCD (dashed lines). 
Using the linear expression the derivatives are calculated at fixed Q 2 > 3 GeV 2 and are 
shown as functions of x in figure [12|. The derivatives show a continuous rise towards low x 
for fixed Q 2 which is well described by the NLO DGLAP QCD to the HI data (section f7T2|). 
The shape of (dF 2 /d\nQ 2 ) x reflects the behaviour of the gluon distribution in the associated 
kinematic range. 



7 QCD Analysis 

In this section the predictions of the DGLAP evolution equations in NLO are confronted 
with the reduced differential cross section measurement. By comparison to the cross section 
data, the strong coupling constant a s , and the shape and normalisation of the gluon and 
quark distributions are determined. This is done using a x 2 minimisation procedure (fit) and 
a suitable choice of parameterisations for the input parton distributions at an initial scale Qq. 

Traditionally, this kind of analysis makes use of both lepton-proton and lepton-deuteron 
data []5%|-|52"[] in order to separate the non- singlet and singlet evolution, and also to deter- 
mine the parton distributions of up and down quarks simultaneously. The present analysis 
introduces a new parameterisation of quark distributions which permits lepton-proton data 
to be analysed alone. Thus the use of deuteron scattering data, which involves bound state 
corrections and their uncertainties, is avoided. 

Two complete analyses are performed, one with HI data only to determine the gluon 
distribution xg(x, Q 2 ) at low x, and a second one in which the HI data are combined with 
data from the BCDMS experiment in order to simultaneously determine the strong coupling 
constant a s (M§) and the gluon distribution. 

The most difficult aspect of these fits is an adequate treatment of systematic errors, which 
lead to strong correlations among the data points. The present analysis uses a sophisticated 
treatment of systematic error correlations allowing their true effect on the extracted quan- 
tities to be estimated by the fit [|52"[]. This procedure is used to identify data regions which 
are strongly affected by correlated systematic errors. The present analysis is therefore based 
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on a minimum number of data sets in regions where their systematic errors are well under- 
stood. Uncertainties due to physics and analysis assumptions are estimated by a systematic 
exploration of the parameter space. 

The Q 2 evolution of parton distributions is the result of processes of radiation from gluons 
which dominate the scaling violations at small x < 0.1, and from quarks which dominate at 
large x. The present HI data allow the gluon contribution to be well determined for fixed a s . 
The strong correlation between the gluon distribution with a s can be much reduced by using 
DIS data at large x and low Q 2 , in addition to the HI data. Thus, the analysis is extended to 
include the precise jip data at large x from the BCDMS collaboration. This allows the gluon 
distribution and a s to be determined simultaneously. 



7.1 Analysis Procedure 

In the quark-parton model, the proton structure function F 2 (x, Q 2 ) is given by a sum of quark 
and anti-quark momentum distribution functions 

F 2 (x, Q 2 )=x^Q 2 g - W, Q 2 ) + q(x, Q% (13) 
g 

where Q q represents the electric charges of quarks. In the present analysis, the sum ex- 
tends over up, down and strange (u, d, s) quarks. The charm and beauty contributions are 



added using NLO QCD calculations [ p3| ] in the on-mass shell renormalisation scheme using 
m c = 1.4 GeV and mj = 4.5 GeV. At low x about 20-30% of the inclusive cross section is 
due to charm production, dominated by the photon-gluon fusion process. 

The present QCD analysis, described in detail in [ffTp, uses a flavour decomposition of 
F 2 (x, Q 2 ) into two independent combinations of parton distribution functions V(x, Q 2 ) and 
A(x, Q 2 ) , according to 

1 11 

F 2 = -xV + — xA. (14) 
3 9 

The x dependences of xg , V and A are parameterised at an initial scale Qfo, and a \ 2 minimi- 
sation determines these distributions and a s . The function V is defined by the valence-quark 
distributions, i.e. 

9 3 

V = -u v - -d v . (15) 



It is thus constrained by the relation 



o 



Vdx = 3, (16) 



which is used in the fit procedure together with the momentum-sum rule. The function 
A contains the sea-quark distribution and a small valence-quark correction. It is given as 

A = u- -(u v -2d v ), (17) 
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and determines the low x behaviour of F 2 (x, Q 2 ) 



These equations hold for a strange contribution s + s = (u + d)/2 and flavour symmetry 
of the sea, u sea = u = d sea = d. This ansatz is generalised in Appendix |A.1| to account 



for the observed small deviations of the strange p4| ] and antiquark [ p5| ] distributions from 



the conventional assumptions about the sea. In the analysis target-mass corrections [q60 are 
applied. 



The analysis is performed in the MS renormalisation scheme using the DGLAP evo- 
lution equations [|]] in NLO [|67|]. Thus the formulae for F 2 given here are modified by 
replacing the sums over parton distributions by sums over integrals of products of coefficient 
functions times parton distributions. 

The strong coupling constant is defined by the solution of the renormalisation group 
equation to order a z s , 

<!<l - -f3 a 2 s -f3 ia l (18) 



where a s = a s /A-K, [i T is the renormalisation scale, and the (3 functions are defined in p7|]. 
The longitudinal structure function is calculated to order a 2 . The analysis uses an x space 
program developed inside the HI collaboration flSBp. This program has been checked in 



detail against different evolution codes [|59|, |69| , |70| ], and very good agreement is found. 



In the fit procedure, a x 2 function is minimised which is defined in Appendix |A.2| . This 



definition takes into account correlations of data points caused by systematic uncertainties. 
It is desirable that the fit results depend neither significantly on the functional form of the 
parameterisations which are used for the input distributions, nor on the input scale Ql at 
which these are defined. Thus, for each fit described below, a grid of nearly 10 3 initial fit 
conditions is considered, with about tenfold variations of each of the parameters Q%, Q 2 min 
and a s . Here Q 2 min denotes the smallest Q 2 value of data included in the fit. 

The quality of the fits is studied in a statistical evaluation of the parameter space for 
all data sets and parton distribution parameterisations considered. As described in Ap- 
pendix |A.3| , this leads to a best choice of 

xq(x) = a g x bq (l - x) Cq [l + d g y/x + e q x), (19) 

for the parameterisations of the functions V, A and xg. The standard value of Ql is 4 GeV 2 , 
but it can be varied over a reasonable range without significantly influencing the result. 



7.2 Fit to HI Data and Determination of the Gluon Distribution 

The measurements of the reduced differential cross section presented here are combined 



with recent data of the HI collaboration [ |T2| ] from the same data taking period, which cover 
the large x range at high Q 2 > 150 GeV 2 . A cut Q 2 < 3000 GeV 2 is applied to eliminate 
the region where electroweak interference effects are important, which involve the structure 
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function xF 3 . Since the HI data have still limited precision at large x, the parameter d g in 
the expression for the gluon distribution (equation \T% is superfluous, see Appendix |A.3| . For 



the fits to the HI data alone it is therefore set to d g = 0. 

The standard fit assumes a fixed a s (M|) = 0. 1 15 and uses all HI data for Q 2 > 3.5 GeV 2 . 
In this range derivatives of F 2 with respect to In Q 2 are measured and found to be described 
by the DGLAP evolution equations (see section |6|). The momentum fraction carried by the 
gluons is 0.43 ± 0.02 (exp) at the input scale Qq = 4 GeV 2 where the error is due to the 
measurement uncertainties. The variation of the gluon momentum fraction with Q 2 is shown 
in figure [13]. The result agrees with previous determinations at Q 2 = 7 GeV 2 []5^,[7l]]. The 
stability of this result has been checked by adding up and also fid data of the BCDMS collab- 
oration. As is shown in figure [13], no significant change is observed. The fit is also repeated 
without using the constraint given by the momentum sum rule. This fit determines the in- 
tegral Jq x(E + g)dx to be 1.016 ± 0.017(ea:p), where £ is the singlet parton distribution 
function, see Appendix |A.3| . This value is found to be nearly independent of the minimum 
Q 2 value of the data included in the analysis. 

The structure function F 2 (x, Q 2 ) is extracted from the reduced cross section data using 
the prediction of the fit for the longitudinal structure function F L (x, Q 2 ) . The result is shown 
in figures [14| and [15|. The data are compared to published jip data of the fixed target muon- 



proton scattering experiments BCDMS and NMC. The solid lines give the result of the QCD 
fit with Q 2 min = 3.5 GeV 2 to the HI data. This fit also describes the fixed target data in 
the non-overlapping regions rather well, except for the data points at x = 0.65 where the 



fit curve is below the BCDMS data. The HI data at this value of x [ |12| ] have a correlated 
systematic uncertainty of 12%, due to the energy scale error for the scattered positron, which 
accomodates the observed difference. 

The x range is restricted at small x by the choice of Q 2 nin - An extension of the analysis 
to low values of Q 2 and x is of interest to study possible deviations from NLO DGLAP 
evolution. The dependence of the fit result on the chosen Q 2 min is studied systematically. 
Figure [T5| shows the HI F 2 data for x < 8 • 10~ 4 together with the fit curves for different 
values of Q 2 min . The fit with Q 2 min = 1.5 GeV 2 describes all the data very well. If Q 2 min is 
raised, the fit curves extrapolated below Q 2 = Q 2 min tend to undershoot the data excluded 
from the fit. The gluon distributions at Q 2 = 5 GeV 2 obtained from these fits are shown in 
figure [H] in the low x range where the gluon distribution is constrained^ They are consistent 
within the estimated uncertainty in the overlapping regions. Extension of this study into 
the region of Q 2 ~ 1 GeV 2 is of interest. It requires precision data in a large range of x. 
For such Q 2 values the gluon distribution xg(x, Q 2 ) , in twist 2 NLO QCD, is observed to 
approximately vanish at low x aquiring a valence-like shape. 

The fit to the HI data determines a s (M|) to be 0.115 with an experimental error of 
±0.005 and an optimum x 2 of 180 for 224 degrees of freedom. This is the first measurement 
of a s with HERA inclusive cross section data alone. The result for the gluon distribution 



1 The measurement of the slope (di^/c* In Q 2 )^ requires at least two data points with different Q 2 above 



Qmin f° r fixed x. Thus the minimum value x m i n at which this slope can be measured depends on Q 2 



It has been shown in [ p5p that the derivative (dFi/dhxQ ) x determines the gluon distribution at a value of 
approximately 2x. The gluon distributions in figure UT\ are therefore shown only down to x ~ 2x ri 



L mill • 
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obtained from the HI measurements is shown in figure [T8|. The innermost error band is due 
to the experimental measurement uncertainty, which for 3 ■ 10~ 4 < x < 0.1 is about 5% and 
decreases to about 3% at Q 2 = 20 GeV 2 . The middle error band illustrates the effect of an 
a s uncertainty of ±0.0017, which is derived in the fit to the combined HI and BCDMS data 
(see below). The outer error band includes in addition the uncertainty of the QCD model as 
described subsequently. For values of x > 0.1 this analysis is not able to reliably determine 
the shape of the gluon distribution because in this range the scaling violations are dominated 
by quark radiation rather than gluon radiation. Yet, the integrated momentum fraction carried 
by gluons at large x is constrained by the momentum sum rule as discussed above. The fit 
determines this contribution to be xg(x, Q 2 )dx = 0.13 ± 0.04 (exp) for Q 2 = 4 GeV 2 . 

The result for the gluon distribution depends on the theoretical framework since xg is 
not an observable. For example, if the massive quark description for charm and beauty 
production is replaced by the massless treatment of heavy quarks, the gluon distribution 
changes as illustrated in figure [T3|. The gluon distribution in the massless fit is about 15% 
lower at small x compared to the standard result. A consistent cross check of this massless 
fit result is obtained with a Mellin n space program [f72|]. 



The analysis of the longitudinal structure function Ft (section 5.2) uses a QCD fit to the 
reduced cross section for y < 0.35. This fit follows exactly the same procedure as that de- 
scribed above. It results in a % 2 of 151 for 180 degrees of freedom and agrees very well with 
the full fit. In particular, the gluon distribution obtained in this fit is nearly indistinguishable 
from xg obtained in the standard fit covering the full y range which is sensitive to Ft ■ Thus, 
xg(x, Q 2 ) appears to be determined by the scaling violations of F 2 (x, Q 2 ) . 



7.3 Simultaneous Determination of a s (M|) and the Gluon Distribution 



The precision of the large x, high Q 2 HI data [ |T2| ] is not sufficient to enable a competitive de- 
termination of a s (M§) and of the gluon distribution simultaneously from the HI data alone. 
The most precise measurement of the DIS inclusive cross section at large x was obtained by 



the BCDMS fip scattering experiment [|T3J (figure [15]). These data are therefore combined 
with the HI measurements. 

In a first step, a fit is performed to the complete data sets. The correlated systematic 
errors of the data are fitted, together with the other parameters. Regions of data are identified 
in which the fit causes large systematic shifts of the data points. For the BCDMS data in the 
range yscDMS = Vn < 0.3 the data points are shifted by more than the quoted systematic 
error. The low y region in this experiment is particularly strongly affected by the energy 
scale uncertainty of the scattered muon, which leads to correlated systematic errors oc 1/y. 
These become large at low y for each of the four data sets at different muon energy. Note 
that the low y data of BCDMS differ from measurements of the ep scattering cross section at 
SLAC [[73|,[7?p in the region of overlap. In this region the BCDMS data accuracy is dominated 
by systematic errors while the SLAC measurement is statistically limited. This also suggests 
the presence of large systematic effects in the low y region of the BCDMS data which were 
studied previously [|75j]. Therefore in all analyses only BCDMS data with y^ > 0.3 are used. 
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The result of the QCD fit to the combined HI and BCDMS data sets, leaving a s as a free 
fit parameter, is shown in figures [14] and [15| as dashed curve. It describes the data very well 
with a x 2 / d °f of 394/451. The fit to the HI and BCDMS data is nearly indistinguishable 
from the fit to the HI data alone, except for the two highest x bins. The parameterisations 
used in the fit to the HI and BCDMS data are given in equation [I9L and the parameters are 
summarised in table |5]. The choice of these particular shapes results from a detailed analysis 
of the behaviour of the x 2 function, similar to the fit to the HI data which is described in 
Appendix [A.3 . 





a 


b 


c 


d 


e 


gluon 


1.10 


-0.247 


17.5 


-4.83 


68.2 


V 


86.3 


1.47 


4.48 


-2.12 


1.60 


A 


0.229 


-0.130 


19.7 


-3.82 


29.8 



Table 5: Parameters of the input distributions xq(x) = a q x bq (1 — x) Cq [1 + d q y/x + e q x] for 
xg(x, Q 2 ) , V(x, Q 2 ) and A(x, Q 2 ) at the initial scale Q 2 = 4 GeV 2 using HI and BCDMS 
data for Q 2 > 3.5 GeV 2 and > 0.3. A fit with d g = e g = yields c g = 6.5, not far from 
the dimensional counting rule expectation jf76[/, yet with a worsened x 2 - 



The a s value obtained in the NLO analysis of the HI and BCDMS proton data is 

a s (M 2 z ) = 0.1150 ± 0.0017 (exp) + Slooos (model). (20) 

In this combined fit both data sets consistently favour a value of a s (M\) ~ 0.115 with 
comparable accuracy. The first error represents the experimental uncertainty of the data 
sets. The second error includes all uncertainties associated with the construction of the QCD 
model for the measured cross section. These are summarised in table ^|. The strong coupling 
constant is defined here by the solution of the renormalisation group equation to order a^. 
In the double logarithmic approximation the value for ol s (M\) is calculated to be lower by 
0.0003. 

The value obtained for cn s (Mf ) is nearly independent of Ql and of the chosen parameter- 
isation for the large set of input distributions considered in Appendix |A.3| . Residual effects 
are included in the estimation of the systematic error on a s . The dependence of a s on Q 2 min 



is shown in figure [20]. No systematic trend is observed. Note that the BCDMS data are lim- 
ited to Q 2 > 7.5 GeV 2 , such that a choice of Q 2 min below this value affects the low i HI 
data only. 

The combination of low x data with high x data constrains the gluon distribution and 
a s . A correlation is observed (figure |2l]) between a s and the parameter b g , which governs 



the shape of the gluon distribution at low x (equation [T^). In a fit to the BCDMS data alone, 
for > 0.3 and using xg = ax b (l — x) c , a value of a s (M|) = 0.111 ± 0.003 (exp) is 
found 0, and b g is positive. A positive value of b g implies that xg(x, Ql) falls as x decreases. 

2 The requirement > 0.3 causes a s (M|) to increase by about 0.004 in the fit to the BCDMS data only, 
and by about 0.002 in the fit to the HI and BCDMS data combined. 
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analysis uncertainty 


+5 a s 


-5 a s 


Q 2 min = 2 GeV 2 

w I l III U 




0.00002 


Q 2 min = 5 GeV 2 

1 1 bill U 


0.00016 




parameterisations 


0.00011 




Ql = 2.5 GeV 2 


0.00023 




Ql = 6 GeV 2 




0.00018 


y e < 0.35 


0.00013 




x < 0.6 


0.00033 




> 0.4 


0.00025 




x > 5 • 10" 4 


0.00051 




uncertainty of u — d 


0.00005 


0.00005 


strange quark contribution e = 


0.00010 




m c + 0.1 GeV 


0.00047 




m c — 0.1 GeV 




0.00044 


m b + 0.2 GeV 


0.00007 




m b - 0.2 GeV 




0.00007 


total uncertainty 


0.00088 


0.00048 



Table 6: Contributions to the error of a s (Mf ) in the analysis of HI ep and BCDMS jip data 
which are due to the selection of data and to the fit assumptions. 



An early a s analysis J60|], in the absence of detailed information about the low x behaviour 
of xg , assumed b g = 0. A positive or zero value of b g , for Ql > 4 GeV 2 , is incompatible 
with analyses of the HERA DIS data at low x. The fit to the BCDMS data, when comple- 
mented with low x HI data, leads to a negative b g parameter, and therefore a larger value of 
a s (M|) is obtained with the BCDMS data than hitherto. The results of these fits to the HI 
and BCDMS data are shown in figure j2"2]a. In the combined fit both data sets give a consistent 
and comparable contribution to the error on a s . This is illustrated in figure [2"2"|b. 

A rather large theoretical uncertainty of the NLO analysis results from the choices of the 
renormalisation scale /i 2 = m r -Q 2 (equation [18]), and of the factorisation scale fij = mj-Q 2 
which leads to scale dependent parton distributions. In the MS scheme both scales are set 
equal to Q 2 , i.e. m r = m/ = 1. In the absence of a clear theoretical prescription, the effect 
of both scales on a s is estimated by varying the scale factors m r and mf between 0.25 and 4. 
The results are summarized in table 0. In agreement with previous studies [ |62| ] it is found that 
the renormalisation scale causes a much larger uncertainty on a a (M|) than the factorisation 
scale. Depending on which set of m r and m/ is chosen, the obtained x 2 differs by several 
units. This suggests that the assumed variation of the scales is too large. These scales, 
however, are not considered to represent physics quantities which may be determined in 
the minimisation procedure. The estimated overall uncertainty of about 0.005 on a s (M|) is 
much larger than the experimental error. It is expected to be significantly reduced when 
next-to-NLO calculations become available [ fTO] , [77] ] . Recently an a s analysis of moments of 
structure functions, measured in charged lepton-nucleon scattering, was presented extending 
toNNLO QCD[^8p. 
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m r = 0.25 


m r = 1 


m T = 4 


m f = 0.25 
m f = 1 
rrif = 4 


-0.0038 
-0.0055 


-0.0001 
+0.0005 


+0.0043 
+0.0047 
+0.0063 



Table 7: Dependence of a s (M z ) on the renormalisation and factorisation scales mf and 
m r , respectively, expressed as the difference of a s (M z ) obtained for scales different from 
one and the central value of a s (M z ) =0.1150. The combination nif = 4 and m r = 0.25 
is abandoned since the splitting function term oc In (m r /m/) 2 becomes negative at low Q 2 
which causes a huge increase of x 2 ■ 



The stability of the fit results is checked further with respect to possible changes in the 
analysis procedure: 

• The value of a„(M|) increases by 0.0005 if the correlation due to systematic errors is 
neglected, i.e. if the correlated systematic error parameters are not part of the minimi- 
sation (see Appendix |A. 2| ) . 



In the present analysis, the relative normalisations of the data sets are left free. The 
change imposed by the fit to the BCDMS data is about —1.5% within a total normali- 
sation uncertainty of 3%. The HI data are moved by about 1% within the experimen- 
tal error of 1.7%. Thus the selected HI and BCDMS data are compatible with each 
other. If the fit is repeated with all normalisations fixed then x 2 increases by 26, and 
a s (M 2 z ) increases by 0.0005. 



• If the BCDMS data is replaced by data of the NMC collaboration [J79[], imposing the 
low Q 2 limit of the BCDMS data, a consistent value of a s (M§) = 0.116 + 0.003 (exp) 
is obtained. 

• If the heavy flavour treatment is changed and a massless, four flavour fit performed, 
a s (M|) is enlarged by +0.0003. 

• The addition of the BCDMS deuteron target data, with y M > 0.3, to the HI and 
BCDMS proton data yields a s (M z ) = 0.1158 ±0.0016 {exp), i.e. a s (Mf) increases 
by 0.0008. In this analysis nuclear corrections [ |80~1 ] are applied, and the conventional 
flavour decomposition into valence and sea quarks is used. 

The gluon distribution from the fit to the HI and the BCDMS proton cross section data 
is shown in figure ^3] for Q 2 = 5 GeV 2 . The inner error band represents the experimental 
uncertainty of the determination of xg for a s fixed. This fit, however, simultaneously deter- 
mines xg(x, Q 2 ) and a s , which leads to a small increase of the experimental error of xg as is 
illustrated by the middle error band. The full error band includes in addition the uncertainties 
connected with the fit ansatz, as listed in table |5| for the determination of a s (M z ) . For the 
low x behaviour of xg(x, Q 2 ) these are dominated by the choice of Q 2 min , as is discussed in 
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section [7~2[ The gluon distribution from the combined fit is shown also for two higher Q 2 
values, 20 and 200 GeV 2 . The DGLAP evolution leads to a gluon distribution which rises 
dramatically at small x with increasing Q 2 (figure The inner solid line illustrates the 
behaviour of xg(x, Q 2 ) , as determined with the HI data alone, which is seen to be in very 
good agreement. The fits with deuteron data or NMC data lead to very similar gluon distribu- 
tions. This analysis determines xg from the scaling violations of F 2 . It is more accurate but 
consistent with determinations of the gluon distribution by the HI experiment in charm [ |8T| ] 
and deep-inelastic dijet [|8~2|] production. 



8 Summary 

A new measurement of the deep-inelastic positron-proton scattering cross section is pre- 
sented for squared four-momentum transfers 1.5 < Q 2 < 150 GeV 2 and Bjorken-x values 
3 • 1CT 5 < x < 0.2 which is more precise than previous measurements in this kinematic 
range. The statistical accuracy of the present inclusive cross section measurement is better 
than 1%, for a large part of the kinematic region. The systematic uncertainty has decreased 
to about 3%, apart from the edges of the covered range. This is due to improved detectors in 
the backward region used for the identification and measurement of the scattered positrons 
at low Q 2 . These are an electromagnetic calorimeter (SPACAL) with very good spatial and 
energy resolutions, a drift chamber (BDC) and a silicon tracker (BST). 

The kinematic range is extended down to y = 0.004 such that the present data over- 
lap kinematically with the measurements of muon-proton scattering experiments. This is 
achieved with higher statistics, improved event vertex reconstruction and calibration of the 
forward parts of the calorimeter. The present data agree with the fxp data in the region of over- 
lap within the accuracy of about 7%. The kinematic range is extended also up to y = 0.82 
using track reconstruction in front of the SPACAL in an extended angular range. 

The cross section measurement is used to determine derivatives with respect to Iny and 
to InQ 2 as functions of x and Q 2 . The partial derivative (<9-F 2 /<91n Q 2 ) x is measured in the 
full x, Q 2 range of this measurement. When considered as a function of x at fixed Q 2 , for 
3 < Q 2 < 40 GeV 2 , it is observed to rise continuously towards low x in agreement with 
QCD. 

The partial derivative of the reduced cross section (da r /d In y) Q 2 is used to extract Fi at 
low Q 2 < 10 GeV 2 . This is complemented by a determination of F L at Q 2 > 10 GeV 2 using 
the difference between the measured reduced cross section, ay, and F 2 calculated from an 
extrapolation of a NLO QCD fit to low y data. Thus the longitudinal structure function F L at 
low x is measured more precisely than hitherto, and in a larger Q 2 range. 

A detailed, systematic analysis is presented of the structure function data using the 
DGLAP evolution equations in NLO. The salient features of this analysis are a new formal- 
ism for charged lepton-proton scattering and a comprehensive study of the influence of model 
parameters and parton distribution parameterisations on the results. These are obtained from 
a minimum number of data sets with special emphasis on the treatment of correlated system- 
atic errors. 
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The scaling violations of F 2 , the behaviour of the derivatives and of the longitudinal 
structure function at low x are found to agree with DGLAP QCD. The present precise low x 
data, when combined with high Q 2 data of HI from the same running period, determine the 
shape of the gluon distribution at small x. The gluon distribution xg(x, Q 2 ) is determined at 
Q 2 = 20 GeV 2 to an experimental accuracy of about 3% in the kinematic range 3 ■ 1CT 4 < 
x < 0.1. 

A simultaneous determination of the gluon distribution and of a s is obtained by com- 
bining the low a; data of HI with /ip scattering data of the BCDMS collaboration at large 
x. A value of the coupling constant a s (M|) = 0.1150 ± 0.0017(ea;p)l o^l(model) is ob- 
tained. The value of a s changes by about 0.005, much more than the experimental error, 
if the renormalisation scale is allowed to vary by a factor of four, and to a lesser extent if 
the factorisation scale is changed by the same amount. This uncertainty is expected to be 
reduced significantly in next-to-NLO perturbation theory. 
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A Details of the QCD Analysis 



A.l Flavour Decomposition of F 2 

The structure function F 2 can be written as 

F 2 = - • xU + - • xD, (21) 
9 9 

with U = u + u and D = d + d + s + s (equation [T3|). A modified projection yields 

F 2 = - ■ xS + - • xA. (22) 
9 3 

The sum S = U + D defines a singlet combination of quark distributions which has a Q 2 
evolution coupled to the gluon momentum distribution xg . The difference A = (2U — D)/2> 
defines a non-singlet distribution which evolves independently of xg . Thus F 2 is defined by 
two independent quark distribution functions. 

In this analysis two specific functions V(x, Q 2 ) and A(x, Q 2 ) are chosen which are re- 
lated to U and D according to 

U — -V + 2A (23) 
3 



and 



D = W + 3A. (24) 



The inverse relations defining V and A are 



V=^(3U- 2D) = 9 -u v - h v + ^u- 3(d + s) (25) 



and 



A=-(2D-U)=d + s-^u--u v + l-d v , (26) 
which for the conventional assumption u = d = 2s lead to the relations presented in the 



introduction of the QCD analysis, see section [7J]. In this approximation the V distribution 



vanishes for small x < 0.01. The behaviour for large x is defined by u v . For small x the 
function A is given by the sea distribution A ~ u. 



Recent measurements of Drell-Yan muon pair production at the Tevatron \ p^\ have es- 
tablished a difference between the u and d distributions. Charged current neutrino-nucleon 
experiments determined the relative amount of strange quarks in the nucleon sea to be 

s + s = (- + e) • (u + d), (27) 
28 



with a recent value of e = —0.08 [|83|]. These results lead to modifications of the simple 
assumptions on the seaf]. They have been accounted for by modifying equation |24] according 
to 



D=-V + kA, 



(28) 



which, using equation |23], results in 



2 fc-1 



(kU - 2D) 



(29) 



and E = V + A • (2 + k). Choosing k = 3 + 2e can be shown to remove the strange 
contribution to the function V yielding 



V = - ■ Y^K 3 + 2t >v ~ 2d v + (5 + 2e)(u - d)}, 



(30) 



which coincides with equation O for e = and u — d. Because the integral 5 = J(u — d)dx 
is finite [], this choice of k allows the counting rule constraint (equation [16|) to be maintained 

as 



/' 

Jo 



Vdx = 3 + 6---^^ = v{e,5). 
o 4 1 + e 



(31) 



If this constraint is released in a fit to the HI data, a value of / Vdx = 2.24 ± 0.13(exp) is 
obtained instead of about 2.5 following from equation ^T]. The modified expression for the 
A function in terms of quark distributions becomes 



.4 



4 1 + e 



[4m - (u v - 2d v ) - 5(u - d) + 2e(u + d)]. 



(32) 



For the naive assumptions e = and u = d this yields the approximate relation [17] and A ~ u 
at low x < 0.1. In the analysis these generalised expressions are used for V, its integral and 
A. 



A.2 Definition of Minimisation Procedure 

The x 2 is computed as 



2 [ a rexp (fj %t x (1 V exp 5 exp J2k^k at ( s k P ))] 2 



es;p dat ^dat,sta ® "dot, unc 

+ E4+EE(^) 2 (33) 



exp exp k 



3 The evolution of s + s in DGLAP QCD is found to yield a linear dependence of e on In Q 2 which is used 
to extrapolate the NuTeV result Q, obtained at 16 GeV 2 , to Q 2 = Qo- 

4 The most accurate measurement of (u— d)dx has been performed by the E866/NuSea Collaboration [ |65| ] 
which obtained a value of -0.118 ± 0.011 at (Q 2 ) = 54 GeV 2 . 



29 



The first two sums run over the data (dat) of the various experiments (exp). 5 exp is the 
relative overall normalisation uncertainty. a dat:Sta and <J Aat ^ unc are the statistical error and 
the uncorrelated systematic error, respectively, corresponding to the datum dat. u exp is the 
number of standard deviations corresponding to the overall normalisation of the experimental 
sample exp. 5 < k lat (s e k xp ) is the relative shift of the datum dat induced by a change by s e k xp 
standard deviations of the k th correlated systematic uncertainty source of the experiment 
exp. 



A.3 Parameterisations 

As explained above, three parton distributions (xg, V and A) are necessary to describe the 
proton structure function F 2 (x, Q 2 ) and its Q 2 evolution. The following general type of 
parameterisation is used 

xq = a q x bq (l - x) Cq [l + d q ^[x + e q x + f q x 2 ] (34) 

for q = g, V and A. An attempt is made to describe these functions with the least number of 
parameters in the brackets of equation R4L All distributions are first calculated using recent 



parameterisations of parton distributions, by GRV98 [|84|], MRS99 [|74|] and CTEQ5 [|85|]. 
The resulting functions are fitted using the expressions of equation ^ in order to obtain 
initial information about how the new linear combinations of parton distributions V and A 
are possibly parameterised best. All global analysis distributions require the presence of dy 
and ey but allow fy to be set to zero. This defines the parameterisation of V which mainly 
is a combination of valence quark distributions, see equation BO. For xg and A, however, 



different parameter combinations are tested in a systematic way using fits to data. 

The choice of a set of parameterisations is guided by the desire for a weak dependence 
of the x 2 function on the initial scale Q%, and by the observed saturation of the x 2 when the 



number of parameters becomes too large. This is demonstrated in figure [24] for the fit to HI 
data alone. The functions without a term oc y/x in the A distribution, see table yield a 
steady decrease of x 2 with Q 2 ,. 

Stability is observed for Ql > 4 GeV 2 for the other parameterisations. Three of them 
have a similar x 2 ■ For the HI fit the parameterisation CP3 is chosen. The functions CP4 and 
CP8 have one more parameter but only one unit of x 2 is gained which points to saturation 
of the parameter list. Although for CP7 a somewhat better x 2 is found, this parameteri- 
sation is not considered. It yields a too large gluon momentum as compared to the other 
fits performed, including those using HI and BCDMS data. These all agree for the gluon 
momentum fraction among each other and also with a previous analysis by the NMC collab- 
oration [|5B|]. In the CP7 fit to the HI data apparently too many parameters are assigned to 
describe the quark distributions at large x since this leads also to a distorted V distribution. 

The choice of parameterisation depends on the data set considered. In a similar study for 
the fit to HI and BCDMS data, the parameterisation CP4 is chosen. Use of parameterisations 
with a high x term (1 + ax 13 ), as introduced by the CTEQ collaboration [|8"5|], worsens the x 2 
by eleven units and has thus not been considered further. 



30 



type 


gluon 


A 


CP1 


1 + ex 


1 + ex 


CP2 


1 + dy/x + ex 


1 + ex 


CP3 


1 + ex 


1 + dy'a; + ex 


CP4 


1 + d^fx + ex 


1 + d-^/x + ex 


CP5 


1 + ex 


1 + ex + fx 2 


CP6 


1 + d^/x + ex 


1 + ex + fx 2 


CP7 


1 + ex 


1 + dy/x + ex + fx 2 


CP8 


1 + d^fx + ex 


1 + dy/x + ex + fx 2 



Table 8: Types of parameterisations of the xg and A distributions at the initial scale Qq. 
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Figure 1 : Distributions of a) the energy, b) the polar angle of the scattered positron, and c) 
Uh for the data sample A taken in 1996/97 (solid points). The histograms show the simulation 
of DIS and the small photoproduction background (shaded), normalised to the luminosity of 
the data. 
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Figure 2: Distributions of a) the energy, b) the polar angle of the scattered positron, and 
c) yh for the low Q 2 data sample B taken in 1997. The histograms represent the simulation 
of DIS and the small photoproduction background (shaded), normalised to the luminosity of 
the data. 
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Figure 3: Distributions illustrating the cross-section measurement at high y (0.46 < y < 
0.82) and low Q 2 (2 < Q 2 < 5 GeV 2 ) for events in the BST acceptance range. DIS event 
distributions of a) the polar angle a) and b) the SPACAL energy of the scattered positron 
candidate, c) SPACAL energy distribution for tagged photoproduction events fulfilling the 
DIS event selection criteria, apart from the E—p z requirement. Solid points: HI data; shaded 
histograms: simulation of photoproduction events; open histograms: added distributions of 
simulated DIS and photoproduction events. 




150 160 5 10 15 10 20 



QJ deg E e '/GeV E e '/GeV 

Figure 4: Distributions illustrating the cross-section measurement at high y (0.46 < y < 
0.89) and large Q 2 (10 < Q 2 < 35 GeV 2 ). a) Polar angle and b) SPACAL energy distribu- 
tions before subtraction of the photoproduction background using the charge measurement 
by the CIC. Solid points: data with positive charge assignment. Shaded histogram: data with 
negative charge assignment. Open histogram: sum of data with negative charge assignment 
and DIS event simulation, normalized to the data luminosity, c) Spectrum of energy mea- 
sured in the electron tagger for DIS candidate events with a linked track of either positive 
charge (solid points) or of negative charge (histogram). 
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Figure 5: Division of the (x, Q 2 ) plane for the measurement of the inclusive DIS cross 
section. At low y the bin size increases as the resolution deteriorates. At large y the data are 
binned in intervals of Q 2 and y in order to account for the y dependent effect of Fl on the 
cross-section and the variation of the systematics with y. The triangular regions inside the 
acceptance do not represent valid analysis bins. 
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Figure 6: Measurement of the reduced DIS scattering cross section (closed points). Triangles 
(squares) represent data from the NMC (BCDMS) muon-proton scattering experiments. The 
solid curves illustrate the cross section obtained in a NLO DGLAP QCD fit to the present 
data at low x, with Q 2 min = 3.5 GeV 2 , and to the HI data at high Q 2 . The dashed curves show 
the extrapolation of this tit towards lower Q 2 . The curves are labelled with the Q 2 value the 
data points belong to and scale factors are coveniently chosen to separate the measurements. 
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Figure 7: Measurement of the derivative da r /dlny at fixed Q 2 . The inner error bars 
represent the statistical errors and the total error bars the statistical and systematic errors, 
added in quadrature. The curves represent the QCD fit result to the HI data, for y < 0.35 
and Q 2 > 3.5 GeV 2 , calculated with different assumptions about Fl . The solid curves use 
the QCD prediction of F L , the dashed (dashed dotted) curves assume F L = F 2 (F L = 0). The 
inner error band is the experimental uncertainty of the fit result (section the outer band 
represents the additional uncertainty due to the fit assumptions. The fit results shown for 
Q 2 < 3.5 GeV 2 (dashed) are obtained by backward extrapolation. 
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Figure 8: Measurement of the reduced DIS scattering cross section (closed points). Triangles 
represent data from the BCDMS muon-proton scattering experiment. The curves represent 
a NLO QCD fit to the HI data alone, using data with y < 0.35 and Q 2 > 3.5 GeV 2 . The 
dashed curves show the F 2 structure function as determined with this fit. The error bands 
represent the experimental and model uncertainty of the QCD fit. 
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Figure 9: The longitudinal structure function F L (x, Q 2 ) for different bins of Q 2 as obtained 
by HI at low x, and by charged lepton-nucleon fixed target experiments at large x. The 
measurements for Q 2 < 10 GeV 2 are determined with the derivative method while the points 
for larger Q 2 are due to the extrapolation method. The error on the data points is the total 
uncertainty of the determination of F L representing the statistical, the systematic and the 
model errors added in quadrature. The inner error bars show the statistical error. The error 
bands are due to the experimental (inner) and model (outer) uncertainty of the calculation of 
F L using the NLO QCD fit to the HI data for y < 0.35 and Q 2 > 3.5 GeV 2 . 



42 



1,5 
1 

0.5 


1.5 
1 

0.5 


1.5 
1 

0.5 



: X = ,0.00008 

— / 

y 

n 


/ 

x = 9<6ooi3 

/ 




_ / 

: X = 0.0002 

/ 

~ 


: X = 0.00032 

— / 

~ 


; X = 0.0005 




: X = 0.0008,'' 

— / 

y 




; X = 0.0013 




: X = 0.002 




: X = 0.0032 




: X = 0.005 




: X = 0.008 




•H1 96-97 
DQCD fit (H1 
P 2 (lnQ 2 ) 



c 
o 



"o 
O 



1 







10 



10 

QV 



Figure 10: Measurement of the proton structure function F 2 (x, Q 2 ) , plotted as functions 
of Q 2 in bins of x (points). The solid lines represent fits to F 2 in bins of x according to a 
polynomial P 2 (x, Q 2 ) = a(x) + b(x) InQ 2 + c(x)(lnQ 2 ) 2 . The dashed lines are obtained 
from the NLO QCD fit to the HI data (section |72|), for Q 2 > 3.5 GeV 2 . The error bands are 
due to the experimental and model uncertainties in the QCD fit. 
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Figure 1 1 : Measurement of the partial derivative (dF 2 /d In Q 2 )x taken at fixed x and plotted 
as functions of Q 2 . The error bars represent the quadratic sum of statistical and systematic 
errors. The straight solid lines are given by the function b(x) + 2c(x) In Q 2 determined in 
fits to F 2 (x, Q 2 ) at fixed x. The dashed lines represent the derivatives as calculated with the 
QCD fit to the HI data. The error bands are due to the experimental and model uncertainties 
in the QCD fit which includes data for Q 2 > 3.5 GeV 2 . 
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Figure 12: The derivative (dF 2 /d\iiQ 2 ) x plotted as functions of x for fixed Q 2 , for the HI 
data (symbols) and the QCD fit to the HI data, for Q 2 > 3.5 GeV 2 (solid lines). The dashed 
curves extrapolate this fit below Q 2 min and outside the range of x. The error bands represent 
the model uncertainty of the QCD analysis. 
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Figure 13: The fraction of the proton momentum carried by gluons as a function of Q 2 , 
obtained in different NLO DGLAP tits. Solid curve: tit to HI data alone; dashed curve: tit 
to HI and BCDMS proton data; dotted curve: tit to HI ep and BCDMS fxp and fxd data. 
The shaded error band represents the experimental uncertainty in the analysis of the HI data 
alone. The solid point is due to a QCD analysis by the NMC collaboration j^Bj. 
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Figure 14: Measurements of the proton structure function F 2 (x, Q 2 ) by the HI and the NMC 
experiments. Solid curves: NLO DGLAP QCD fit to the HI cross section data. Overlayed as 
dashed curves are the results of the QCD fit to the HI ep and BCDMS fip data, for > 0.3, 
which are indistinguishable from those of the pure HI fit. Dotted curves: tit extrapolations 
at fixed x into the region below Q 2 = 3.5 GeV 2 . 
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Figure 15: Measurements of the proton structure function F 2 (x, Q 2 ) by the HI experiment 
and by fixed target muon-proton scattering experiments. The error on the data points is 
the total mesurement uncertainty. The inner error bars represent the statistical error. Solid 
curves: fit to the HI cross section data. Dashed curves: fit to the HI ep and BCDMS fip data, 
fory^ > 0.3. Dotted curves: extrapolations to data not used in the fit. 
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Figure 16: Effect of the Q 2 min cut on the structure function F 2 in the DGLAP QCD fit to the 
HI data (points). The curves represent fits with different minimum Q 2 values. The analysis 
uses Q 2 min = 3.5 GeV 2 as default. 
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Figure 17: Eifecf of the Q 2 min cut, applied in the DGLAP QCD fit to the HI data, on 
the gluon distribution at Q 2 = 5 GeV 2 . The distributions are shown down to x values 
corresponding to twice the minimum x values of the data which allow a Q 2 slope to be 
measured. 
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Figure 18: The gluon distribution xg(x, Q 2 ) at Q 2 = 5 GeV 2 , determined in the NLO 
DGLAP QCD fit to the HI data. Inner error band: experimental uncertainty; middle error 
band: effect of the experimental error and of the a s {Ml) uncertainty of ±0.0017; outer error 
band: effect of experimental, a s and model uncertainties. 
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Figure 19: Comparison of gluon distributions obtained in NLO DGLAP QCD fits to the HI 
data, using different prescriptions: solid curve: standard fit using the massive heavy flavour 
scheme; dashed curve: fit in the massless scheme; dotted curve: fit in the massless scheme 
using a Mellin n space program. 
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Figure 20: Dependence of a s {M\) obtained in fits to the HI and BCDMS data on the min- 
imum Q 2 value used. The error bars denote the experimental uncertainty of a s (M|) . Note 
that the BCDMS data have an intrinsic Q 2 min of 7.5 GeV 2 and are limited in this analysis to 
Un > 0.3 (see text). An increase of Q 2 min implies that the minimum x rises correspondingly, 
i.e. from x = 3.2 • 1(T 5 at Q 2 min = 1.5 GeV 2 to x = 8 • 10" 4 at Q 2 . = 12 GeV 2 . 
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Figure 21 : The correlation of the gluon distribution parameter b g with a s in the fit to HI and 
BCDMS data. The parameter b g governs the low x behaviour of xg oc x ha . 
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Figure 22: Determination of the strong coupling constant a s (M|) in NLO DGLAP QCD. a) 
total x 2 for fits to the HI ep and BCDMS jip data (y^ > 0.3) separately and for the fit using 
data of the two experiments combined; b) partial x 2 contributions of the HI and BCDMS 
proton data in the fit to determine a s using both experiments. 
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Figure 23: Gluon distribution resulting from the NLO DGLAP QCD fit to HI ep and 
BCDMS fip cross section data in the massive heavy flavour scheme. The innermost error 
bands represent the experimental error for fixed a s (M|) =0.1150. The middle error bands 
include in addition the contribution due to the simultaneous fit of a s . The outer error bands 
also include the uncertainties related to the QCD model and data range. The solid lines inside 
the error band represent the gluon distribution obtained in the fit to the HI data alone. 
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Figure 24: Dependence of x 2 on the initial scale parameter Ql for different parameterisa- 
tions of the parton distributions xg and A (appendix [A.3| , table in the NLO QCD fit to the 
HI data. 
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Table 9: The HI measurement of the reduced deep-inelastic cross section with data taken 
in a dedicated low Q 2 trigger run in 1997. Fory < 0.6 the structure function F 2 is extracted 
using the quoted values of R, derived from a QCD fit to the HI cross section data. Frac- 
tional cross section errors are quoted in %. The total error (5 to t) is the quadratic sum of the 
uncorrected (5 unc ), the correlated (5 cor ) and the experimental statistical error (5 sta ). 
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Table 10: The HI measurement of the reduced deep-inelastic cross section with data taken 
in a dedicated low Q 2 trigger run in 1997. Fory < 0.6 the structure function F 2 is extracted 
using the quoted values of R, derived from a QCD fit to the HI cross section data. Fractional 
cross section errors are quoted in %. 
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0.575 


0.244 


0.575 


2.7 


0.7 


1.4 


2.2 


15.0 


0.0320 


0.005 


0.453 


0.211 


0.453 


6.4 


1.0 


3.3 


5.5 


20.0 


0.000268 


0.825 


1.313 






5.2 


3.2 


3.6 


1.8 


20.0 


0.000328 


0.675 


1.383 






2.7 


1.0 


2.0 


1.5 


20.0 


0.000500 


0.443 


1.285 


0.246 


1.324 


2.0 


0.6 


1.6 


0.9 


20.0 


0.000800 


0.277 


1.178 


0.241 


1.190 


1.8 


0.6 


1.4 


1.0 


20.0 


0.00130 


0.170 


1.059 


0.237 


1.062 


1.8 


0.6 


1.4 


1.0 


20.0 


0.00200 


0.111 


0.939 


0.235 


0.940 


2.8 


0.7 


1.4 


2.3 


20.0 


0.00320 


0.069 


0.819 


0.234 


0.819 


2.3 


0.7 


1.5 


1.7 


20.0 


0.00500 


0.044 


0.747 


0.235 


0.747 


2.2 


0.8 


1.5 


1.3 


20.0 


0.0100 


0.022 


0.610 


0.238 


0.610 


2.4 


0.7 


1.4 


1.8 


20.0 


0.0320 


0.007 


0.455 


0.198 


0.455 


6.8 


1.0 


3.3 


5.9 


25.0 


0.000335 


0.825 


1.379 






5.9 


4.1 


3.9 


1.8 


25.0 


0.000410 


0.675 


1.371 






2.6 


1.2 


2.0 


1.2 


25.0 


0.000500 


0.553 


1.345 


0.248 


1.417 


2.4 


1.0 


1.8 


1.2 


25.0 


0.000800 


0.346 


1.242 


0.243 


1.263 


1.9 


0.7 


1.6 


0.9 


25.0 


0.00130 


0.213 


1.091 


0.238 


1.097 


1.8 


0.7 


1.4 


0.9 


25.0 


0.00200 


0.138 


0.985 


0.236 


0.987 


2.9 


0.8 


1.4 


2.4 


25.0 


0.00320 


0.086 


0.879 


0.234 


0.880 


2.8 


0.8 


1.5 


2.2 


25.0 


0.00500 


0.055 


0.754 


0.234 


0.754 


2.4 


0.9 


1.5 


1.6 


25.0 


0.00800 


0.034 


0.663 


0.234 


0.663 


2.5 


0.9 


1.5 


1.8 


25.0 


0.0158 


0.018 


0.547 


0.226 


0.547 


3.7 


0.9 


1.5 


3.3 


25.0 


0.0500 


0.005 


0.447 


0.148 


0.447 


7.5 


1.3 


3.3 


6.6 


35.0 


0.000574 


0.675 


1.473 






2.7 


1.4 


2.0 


1.2 


35.0 


0.000800 


0.484 


1.354 


0.244 


1.405 


2.2 


0.9 


1.7 


1.1 


35.0 


0.00130 


0.298 


1.181 


0.239 


1.195 


1.8 


0.8 


1.4 


0.8 


35.0 


0.00200 


0.194 


1.031 


0.235 


1.035 


1.8 


0.8 


1.4 


0.8 



Table 11: The HI measurement of the reduced deep-inelastic cross section with data taken 
in 1996/97. For y < 0.6 the structure function F 2 is extracted using the quoted values of 
R, derived from a QCD tit to the HI cross section data. Fractional cross section errors are 
quoted in %. 
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QVGeV 2 


X 


V 




R 


F 2 


Stot 


gta 






35.0 


0.00320 


0.121 


0.935 


0.233 


0.936 


3.1 


0.9 


1.5 


2.5 


35.0 


0.00500 


0.077 


0.821 


0.231 


0.821 


2.7 


0.9 


1.5 


2.0 


35.0 


0.00800 


0.048 


0.719 


0.228 


0.719 


2.4 


1.0 


1.5 


1.6 


35.0 


0.0130 


0.030 


0.625 


0.222 


0.625 


3.2 


1.1 


1.6 


2.6 


35.0 


0.0251 


0.015 


0.524 


0.195 


0.524 


4.1 


1.1 


1.6 


3.6 


35.0 


0.0800 


0.005 


0.413 


0.095 


0.413 


9.2 


1.8 


3.5 


8.3 


45.0 


0.00130 


0.383 


1.282 


0.238 


1.309 


1.9 


0.9 


1.7 


0.3 


45.0 


0.00200 


0.249 


1.107 


0.234 


1.115 


1.8 


0.9 


1.4 


0.6 


45.0 


0.00320 


0.156 


0.979 


0.231 


0.982 


1.8 


0.9 


1.4 


0.7 


45.0 


0.00500 


0.099 


0.872 


0.228 


0.873 


2.8 


1.1 


1.5 


2.1 


45.0 


0.00800 


0.062 


0.743 


0.224 


0.743 


2.5 


1.1 


1.5 


1.6 


45.0 


0.0130 


0.038 


0.649 


0.215 


0.649 


2.8 


1.3 


1.6 


2.0 


45.0 


0.0251 


0.020 


0.525 


0.187 


0.525 


4.3 


1.2 


1.6 


3.8 


45.0 


0.0800 


0.006 


0.396 


0.091 


0.396 


7.6 


2.1 


3.5 


6.4 


60.0 


0.00200 


0.332 


1.245 


0.231 


1.263 


2.1 


1.0 


1.7 


0.7 


60.0 


0.00320 


0.208 


1.052 


0.227 


1.057 


1.9 


1.1 


1.4 


0.7 


60.0 


0.00500 


0.133 


0.900 


0.223 


0.902 


3.1 


1.2 


1.6 


2.3 


60.0 


0.00800 


0.083 


0.803 


0.218 


0.804 


2.8 


1.3 


1.6 


1.9 


60.0 


0.0130 


0.051 


0.683 


0.208 


0.683 


2.9 


1.4 


1.7 


2.0 


60.0 


0.0200 


0.033 


0.597 


0.192 


0.597 


4.1 


1.7 


1.8 


3.2 


60.0 


0.0398 


0.017 


0.506 


0.145 


0.506 


4.7 


1.8 


1.9 


3.9 


60.0 


0.130 


0.005 


0.360 


0.057 


0.360 


9.4 


3.0 


3.8 


8.1 


90.0 


0.00320 


0.311 


1.107 


0.221 


1.120 


2.6 


1.4 


1.8 


1.2 


90.0 


0.00500 


0.199 


0.999 


0.216 


1.004 


2.2 


1.3 


1.5 


0.9 


90.0 


0.00800 


0.124 


0.845 


0.209 


0.846 


3.3 


1.5 


1.7 


2.5 


90.0 


0.0130 


0.076 


0.728 


0.197 


0.728 


2.8 


1.6 


1.8 


1.4 


90.0 


0.0200 


0.050 


0.618 


0.180 


0.618 


3.9 


1.8 


1.9 


2.8 


90.0 


0.0398 


0.025 


0.506 


0.135 


0.506 


3.8 


2.0 


1.9 


2.7 


90.0 


0.130 


0.008 


0.339 


0.053 


0.339 


4.8 


3.3 


2.0 


2.9 


120.0 


0.00500 


0.266 


1.011 


0.210 


1.019 


3.7 


2.4 


2.1 


1.9 


120.0 


0.00800 


0.166 


0.839 


0.202 


0.841 


3.0 


1.9 


1.8 


1.4 


120.0 


0.0130 


0.102 


0.744 


0.190 


0.745 


4.6 


2.1 


2.0 


3.6 


120.0 


0.0200 


0.066 


0.604 


0.173 


0.605 


4.0 


2.3 


2.1 


2.5 


120.0 


0.0320 


0.041 


0.558 


0.145 


0.558 


5.8 


2.8 


2.5 


4.5 


120.0 


0.0631 


0.021 


0.462 


0.094 


0.462 


4.9 


3.0 


2.7 


2.9 


120.0 


0.200 


0.007 


0.312 


0.037 


0.312 


10.6 


4.8 


4.8 


8.2 


150.0 


0.0200 


0.083 


0.709 


0.167 


0.709 


8.9 


4.4 


3.8 


6.6 


150.0 


0.0320 


0.052 


0.550 


0.140 


0.550 


8.0 


5.1 


4.1 


4.7 


150.0 


0.0631 


0.026 


0.418 


0.090 


0.418 


9.0 


5.7 


4.4 


5.4 


150.0 


0.200 


0.008 


0.296 


0.036 


0.296 


12.9 


8.4 


7.1 


6.7 



Table 12: The HI measurement of the reduced deep-inelastic cross section with data taken 
in 1996/97. For y < 0.6 the structure function F 2 is extracted using the quoted values of 
R, derived from a QCD tit to the HI cross section data. Fractional cross section errors are 
quoted in %. 
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Q 2 /GeV 2 


X 


V 


(da r /d\nv)n2 




^■sys 




1.5 


0.000039 


0.425 


0.104 


0.039 


0.080 


0.089 


1.5 


0.000062 


0.269 


0.088 


0.070 


0.114 


0.134 


2.2 


0.000045 


0.538 


0.051 


0.017 


0.061 


0.064 


2.2 


0.000075 


0.325 


0.139 


0.015 


0.025 


0.029 


2.2 


0.000122 


0.200 


0.126 


0.015 


0.024 


0.028 


2.2 


0.000194 


0.125 


0.118 


0.018 


0.031 


0.036 


2.2 


0.000291 


0.083 


0.070 


0.018 


0.026 


0.032 


2.2 


0.000454 


0.054 


0.082 


0.014 


0.026 


0.030 


2.2 


0.000748 


0.032 


0.086 


0.018 


0.033 


0.038 


2.2 


0.00122 


0.020 


0.081 


0.016 


0.039 


0.042 


2.2 


0.00224 


0.011 


0.100 


0.007 


0.023 


0.024 


4.2 


0.000086 


0.538 


0.002 


0.031 


0.057 


0.065 


4.2 


0.000143 


0.325 


0.122 


0.020 


0.033 


0.039 


4.2 


0.000232 


0.200 


0.148 


0.017 


0.023 


0.029 


4.2 


0.000371 


0.125 


0.170 


0.020 


0.027 


0.034 


4.2 


0.000556 


0.083 


0.140 


0.019 


0.019 


0.027 


4.2 


0.000867 


0.054 


0.118 


0.014 


0.018 


0.023 


4.2 


0.00143 


0.032 


0.095 


0.014 


0.016 


0.022 


4.2 


0.00232 


0.020 


0.114 


0.011 


0.013 


0.017 


4.2 


0.00488 


0.009 


0.106 


0.004 


0.013 


0.013 


7.5 


0.000154 


0.538 


-0.004 


0.034 


0.049 


0.059 


7.5 


0.000255 


0.325 


0.218 


0.028 


0.024 


0.037 


7.5 


0.000414 


0.200 


0.193 


0.024 


0.031 


0.039 


7.5 


0.000663 


0.125 


0.182 


0.028 


0.027 


0.039 


7.5 


0.000992 


0.083 


0.203 


0.026 


0.021 


0.034 


7.5 


0.00155 


0.054 


0.127 


0.019 


0.018 


0.026 


7.5 


0.00255 


0.032 


0.189 


0.024 


0.021 


0.032 


7.5 


0.00414 


0.020 


0.079 


0.019 


0.021 


0.029 


7.5 


0.00764 


0.011 


0.149 


0.008 


0.016 


0.018 


7.5 


0.0155 


0.005 


0.081 


0.013 


0.012 


0.018 


13.5 


0.000199 


0.748 


-0.399 


0.167 


0.191 


0.253 


13.5 


0.000277 


0.538 


0.151 


0.018 


0.051 


0.054 


13.5 


0.000459 


0.325 


0.182 


0.015 


0.024 


0.028 


13.5 


0.000746 


0.200 


0.299 


0.016 


0.054 


0.056 


13.5 


0.00119 


0.125 


0.295 


0.037 


0.043 


0.057 


13.5 


0.00179 


0.083 


0.148 


0.037 


0.035 


0.051 


13.5 


0.00279 


0.054 


0.173 


0.012 


0.014 


0.018 


13.5 


0.00459 


0.032 


0.203 


0.019 


0.017 


0.025 


13.5 


0.00746 


0.020 


0.145 


0.018 


0.027 


0.032 


13.5 


0.0137 


0.011 


0.135 


0.010 


0.015 


0.018 


13.5 


0.0279 


0.005 


0.045 


0.015 


0.035 


0.038 



Table 13: The HI measurement of the cross section derivative (da r /d\ny)Q2 = 
— (da r /d\nx)Q2 calculated at fixed Q 2 . The errors are given in absolute values. The data 
below 13.5 GeV 2 belong to the special low Q 2 run taken in 1997. The data at larger Q 2 were 
taken in 1996/97. 62 
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0.048 
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0.352 
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0.071 


135.0 


0.0119 


0.125 


0.176 
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0.292 
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0.041 


0.064 
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Table 14: The HI measurement of the cross section derivative (da r /d\ny)Q2 = 
— (da r /d\nx)Q2 calculated at fixed Q 2 . A sta denotes the uncertainty due to the data statis- 
tics. The uncertainties are given in absolute values. 
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F L (x,Q 2 ) 
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2.2 


0.000045 


0.538 


0.100 


0.030 


0.107 


0.025 


0.114 


4.2 


0.000086 


0.538 


0.273 


0.055 
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0.027 


0.118 


7.5 


0.000154 


0.538 
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0.058 


0.088 


0.039 


0.112 


12.0 
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0.825 


0.429 
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0.045 


0.130 


12.0 
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0.210 
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0.036 


0.118 
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0.035 
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20.0 
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0.039 
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0.029 


0.128 
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0.484 


0.239 


0.064 


0.213 


0.047 


0.228 



Table 15: The HI determination of the longitudinal structure function Fl(x,Q 2 ) . The 
errors are given in absolute, 5 sta representing the experimental statistics. The systematic 
errors consider all contributions from correlated and uncorrelated systematic error sources. 
5 met is due to the uncertainties connected with the representation of F 2 in the derivative 
method, for Q 2 < 10 GeV 2 , and in the QCD extrapolation method, for Q 2 > 10 GeV 2 . 
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0.00050 


0.495 
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0.158 
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0.103 
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10.1 
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22.4 
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0.050 


0.052 


0.072 


29.6 
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0.409 


0.423 


0.045 


0.044 


0.062 



Table 16: The HI Measurement of the derivative (dF 2 /d\nQ 2 ) x . A sta denotes the uncer- 
tainty due to the data statistics. The uncertainties are given in absolute values. 
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Table 17: The HI measurement of the derivative (dF 2 /d\nQ 2 ) x . A sta denotes the uncer- 
tainty due to the data statistics. The uncertainties are given in absolute values. 
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Table 18: The HI measurement of the derivative dF 2 /dlnQ 2 . A sta denotes the uncertainty 
due to the data statistics. The uncertainties are given in absolute values. 
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